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Abstract 



In this paper, the duality between the Stieltjes electrostatic interpretation for zeros 
of orthogonal polynomials and the quantum Hamilton Jacobi formalism is shown. From 
this duality, the bound state problem mimics as n unit moving imaginary charges 
ih, which are placed in between the two fixed imaginary charges arising due to the 
classical turning points of the potential. The interaction potential between n unit 
moving imaginary charges ih is given by logarithm of the wave function. For an exactly 
solvable potential, this system attains stable equilibrium position at the zeros of the 

m ■ 

orthogonal polynomials depending upon the interval of the classical turning points. 
keywords : Orthogonal polynomials, quantum Hamilton Jacobi and zeros of orthogonal 



polynomials. 

1 Introduction 

Stieltjes [TJ |2] considered the following problem with n moving unit charges, interacting 
through a logarithmic potential, are placed between two fixed charges p and q at —1 and 
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1 respectively on a real line. He then proved that the system attains a stable equilibrium 
position at the zeros of the Jacobi polynomial P^^\x). Proof is given in Szego's book 
(section 6.7) [3]. If, the interval is changed on the real line, for the fixed charges, then the the 
system attains stable equilibrium position at the zeros of the orthogonal polynomial with the 
respective intervals. For example, in the interval [0; oo) one gets the Laguerre polynomials 
L^\x) and for the the interval (— oo; oo) one gets the Hermite polynomials polynomials 
H n (x). This model has been extended to the zeros of general orthogonal polynomials in the 
ref @]. 

The Quantum Hamilton Jacobi (QHJ) formalism, was formulated for the bound state 
problems by Leacock and Padgett [5j [6] and later on was successfully applied to several 
exactly solvable models (ESM) [7J [H [9j [TOj E] in one dimension, the quasi - exactly solvable 
(QES) models [12], the periodic potentials p2] and the PT symmetric potentials [2] in 
quantum mechanics. In QHJ the central role is played by the quantum momentum function 
(QMF). This function, in general, contains fixed poles that arises due to the classical turning 
points of the potential. In general, for most of the potentials in quantum mechanics there 
will be only two fixed poles, and n moving poles arise due the zeroes of wave function. Thus, 
one can immediately see the connection between the two scenarios presented above. The 
fixed poles of the potential are like the the two fixed charges and n moving poles on the real 
line are like n moving charges. In this letter the duality between the two models in shown 
in section 2. The brief introduction of these models will follow below. 

1.1 Electrostatic Model 

Stieltjes considered the interaction forces for the n moving unit charges arising from a loga- 
rithmic potential which are in between the to fixed charges p and q at — 1 and 1 respectively 
on a real line as 



L = -LogD n (xi,x 2 ...x n ) + p^Log{ — 
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Then, he proved in ref [HI2] that the expression ([T]) becomes a minimum when {x\, x 2 , ■ ■ ■ , x n ) 
are the zeros of the Jacobi polynomial. For completeness the proof is given as in Szego's 
book (section 6.7) [3j. 

The minimum for all the Xi are distinct and different from ±1 and it occur for dL/dx k = 
0(1 < k < n) so that one gets the system of equations 
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By introducing the polynomial 



f n (x) = {x- xi)(x - x 2 ) ■ ■ ■ (x - x n ), 



then one has 
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and using the L'Hospital's rule 
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then this becomes 
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This means that the polynomial 

(1 - x 2 )f^(x) + 2[q-p-(p + q)x\f n {x) = (9) 

vanishes at the points Xk and since this polynomial is of degree n it must be a multiple of 
f n {x)- The factor is easily obtained by equating the coefficient of x n and we have 

(1 - x 2 )f^x) + 2[q-p-(p + q)x}f n {x) = n[n + 2(p + q) - l]f n (x) (10) 

which is the differential equation for the Jacobi polynomial P^ 2p ~ 1 ' 2q ~ 1 \x) . The zeros of the 
Laguerre and the Hermite polynomials admit the same interpretation. 

1.2 Quantum Hamilton Jacobi 

In this section, a brief review of Quantum Hamilton Jacobi formalism is presented below. 
For details see the references [TTj, [9]. The Schrodinger equation is given by, 

h 2 

- 7T- V 2 ^(x, y, z) + V(x, y, z)ip(x, y, z) = Eip(x, y, z). (11) 
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One defines a function S analogous to the classical characteristic function by the relation 

/iS\ 
^(x,y,z) = exp(—J (12) 

which, when substituted in (fiTj) . gives 

(VS) 2 - ihV.(VS) = 2m(E - V(x, y, z)). (13) 

the quantum momentum function p is defined in terms of the function 5" as 

p = VS. (14) 

Substituting ffT4l in ffT3|) gives the QHJ equation for p as 

(p) 2 -ihV.p = 2m(E-V(x,y,z)) (15) 



and from f ill I) and (|14p . one can see that p is the the logarithmic derivative of ip(x, y, z) i. e. 

p=—ihVlmp(x,y,z) (16) 

The above discussion of the QHJ formalism is done in three dimensions the same equation 
in one dimension takes the following form 

p 2 -ih-?- = 2m(E-V(x)), (17) 

ax 

which is also known as the Riccati equation. In one dimension the eq (TTB1) take the form 

p = —ih—lnip(x). (18) 

ax 

It is shown by Leacock and Padgett [5j [6] that the action angle variable gives rise to exact 
quantization condition 

J(E) = — I pdx = nh. (19) 

2 Duality 

By considering the form of the wave function, in the equation (1T81 . to be ip = Y\^ =1 (x — Xi). 
Then, in the quantum momentum function this corresponds to n zeros on the real line. These 
zeros are also called the moving poles in the language of QHJ. Choosing, an exactly solvable 
potential V(x), with two fixed poles as the classical turning points, substituting in equation 
(IT7|) . Then, for bound states the following feature always arises in QHJ that the n moving 
poles lie in between the two fixed poles and the solutions are the orthogonal polynomials for 
the exactly solvable potential V(x). The examples are the Harmonic oscillator, the Coulomb 
potential , the Scarf potential etc 0, [TTj . 

Thus, the connection between the QHJ and the Stieltjes electro static model can be seen. 
The fixed poles of the potential are like the two fixed charges and the n moving poles of the 
real line are like n moving charges. In the electrostatic model the moving charges interact 



with the logarithmic potential and in the QHJ the logarithmic potential arises from the wave 
function. As the quantum momentum function is log derivative of the wave function. 

Starting with the QMF the duality between the two models is established The fact that 
only the residues of the QMF are required for finding the eigenvalues is studied in ref [7J |8] . 
The formalism for effectively obtaining both the eigenfunctions and the eigenvalues from the 
singularity structure of the quantum momentum function is given in ref |9] . The quantum 
momentum function assume that there are no other singular points of p in the complex plane. 
Then the quantum momentum function is given by [3 El El QUI III 02] 

n -i 
p=J2^— + Q(x), (20) 

fc=i x Xk 
here the moving poles are simple poles with residue — %h (we take here h = m = 1) [91 [TT] 

and Q(x) is the residues of fixed poles arising due to the exactly solvable potentials. This 

equation resembles the equation ([3]) except that it is the minimum of the potential. Thus, the 

quantum momentum function can interpret as system of equations arising for the logarithmic 

derivative of wave function and fixed poles arising from the classical turning points. By 

asking the following question, when does this system come to stable equilibrium ? From the 

above discussion it is clear that answer can be obtained using Stieltjes Electrostatic model. 

It can be shown that the same wave function can be obtained from both the models. Thus, 

their exist a duality between the Stieltjes electrostatic interpretation for zeros of orthogonal 

polynomials and the quantum Hamilton Jacobi formalism. 

By rewriting the equation (1201) and solving it as a system of equations 



n -i 

P=J2^— +iQ(x) = 0, (21) 

fe=i x Xk 



and applying the L'Hospital's rule 
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+ iQ(x) = 0. (23) 



By demanding the points x k to vanish, the solution to equation (1251) will be in terms of 
the orthogonal polynomial for an exactly solvable potentials. The interval is fixed by the 
fixed poles of the potential. It is well known that the classical orthogonal polynomials arise 
as solutions to the bound states problems. Thus, the classical orthogonal polynomials are 
classified into three different categories depending upon the range of the polynomials. The 
polynomials in the intervals (—00; 00) are the Hermite polynomials, in the intervals [0; 00) 
are the Laguerre polynomials and in the intervals [— 1; 1] are the Jacobi polynomials. Their 
singularity structure is as follows Q{x) = x , Q{x) — | + C, and Q(x) = — ^j — jtj for the 
Hermite, the Laguerre and the Jacobi polynomials respectively. Hence, the differential equa- 
tion can be obtained by examining at the singularity structure of the quantum momentum 
function. This can be seen by rewriting the eq ( )23|) as 

ft{x)+2iQ(x)&(x) = 0. (24) 

The function Q(x) which has the information of fixed pole singularity structure appears 
as the coefficient of f' n {xk)- By examining the differential equations of the Hermite, the 
Laguerre and the Jacobi polynomials the coefficients of Q(x) are fixed. 
Let f(x) = L™(x) satisfy the Laguerre differential equation 

x-^ffr) + (m + l- x)£f(x) + Xf(x) = 0, (25) 

where A is an integer. By examining the first two terms of the differential equations ( I24p 
and (125]) one gets 

2iQ{x) = ^±3. _ 1 (26) 

x 

the singularity structure for the Laguerre is 

Q( x ) = - + C (27) 

x 

thus one gets b = —iim + 1) and C = i. Similarly for the Jacobi differential equation 

(1 - x 2 )f:(x k ) + 2[p-q-(p + q)x\f n {x k ) + n[n + 2(p + q) - l]f n (x) = (28) 



again comparing the first two terms 



the singularity structure for the Jacobi is 

^Q(x) = — - -±- (30) 

x — 1 X + 1 

thus one has p = —ia and q = —ib. Similar analysis can be done for the Hermite polynomials. 

The values of m, p and q has to be determined as these are not points like in the 

electrostatic model. The method adopted by QHJ, search for the polynomial solutions leads 

to quantization, are used to calculate these values. By writing the quantum momentum 

function as 

n fir \ 

p = E*7^ + <2(*) ( 31 ) 

and substituting in ([17]) then one gets 

ft{x k ) + 2iQ[x)f' n [x k ) + \Q\x) - iQ\x) -E + V(x)]f(x) = 0. (32) 

by demanding [Q 2 (x)— iQ'(x) — E + V(x)] to be constant, the values of the residues appearing 
for fixed poles are obtained, in the process the system is quantized. Thus by comparing the 
equation (J231) and (1321 it can seen that the singularity structure of iQ(x) determines the 
differential equation. Therefore, the same wave function is obtained from both the methods. 

3 Discussion 

From the previous discussion, it is clear that the two models are dual to each other. Therefore, 
From this duality, the bound state problem mimics as n unit moving imaginary charges ih, 
which are placed in between the two fixed imaginary charges arising due to the classical 
turning points of the potential. The interaction potential between n unit moving imaginary 
charges ih is given by logarithm of the wave function. For an exactly solvable potential, 
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this system attains stable equilibrium position at the zeros of the orthogonal polynomials 
depending upon the interval of the classical turning points. Once charges arise in any model 
they satisfy the continuity equation of the form 

^P + V-J = 0. (33) 

Since, the equation ( |24p and (132 j) are nothing but the different form of the Schroedinger 
equation. Therefore their exist a continuity equation of this form for these imaginary with 
P = Jv^^PdV is probability density function and J = h[ip*(Vip) — ^(V^>*)] is probability 
current density function. Hence, the conservation of probability leading to conservation of 
imaginary charge and probability current leads to current density for imaginary charge. In 
this model p is the amount of imaginary charge and J is the current density for imaginary 
charge. Thus, this model is consistent with quantum mechanics. 

4 Conclusion 

In this paper, the two different models, one the Stieltjes electrostatic model and the other one 
Quantum Hamilton Jacobi formalism are examined. It is shown that the two models are dual 
to each other except that one is a classical model and another is a quantum model. One new 
feature comes out of this study is that the wave function can be obtained from the quantum 
momentum function itself, one need not solve the quantum Hamilton Jacobi equation. From 
Stieltjes electrostatic model gives nice insights to the methodology of quantum Hamilton 
Jacobi formalism. It is interesting to note that the Stieltjes electrostatic model existed 
almost 30 years before quantum mechanics came into existence. 
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